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I  INTRODUCTION 


The  intense  electromagnetic  fields  obtp.inable  in  laser  beams  can 

induce  multiphoton  transitions  in  atoms  and  molecules  with  significant 

* 

probabilities.  A  variety  of  such  transitions  have  been  observed. 

Also,  a  number  of  theoretical  studies  of  these  processes  have  been 
published.  The  simplest  of  these  are  essentially  dimensional  analyses3 
and  are  very  crude,  but  do  provide  useful  scaling  relations  in  the  weak 
field  limit. 

The  most  straightforward  of  the  calculations  of  multiphoton  tran¬ 
sitions  use  perturbation  theory.4  These  calculations  are,  unfortunately, 
of  great  numerical  complexity  and  are  consequently  ordinarily  carried  out 
with  such  gross  approximations  that  the  results  are  probably  not  much 
more  accurate  than  are  those  of  the  above-mentioned  dimensional  analyses. 
The  perturbation  calculations,  of  course,  fail  to  describe  transitions 
in  very  strong  fields. 

A  somewhat  intuitive  and  partly  phenomenological  calculation  by 
Keld”sh6  takes  account  of  distortion  of  the  atom  or  molecule  by  the 
intense  fields  and  also  of  transitions  taking  place  through  a  nearly 
resonant  intermediate  state.  Keldysh's  theory  is  not  restricted  to  the 


* 

Several  reviews  of  multiphoton  processes  are  available  in  the 
literature.1^  See  also  R.  A.  Fox  ct  al.2 

^  References  are  listed  at  the  end  of  the  report. 
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weak  field  limit,  as  is  perturlation  theory,  and  at  the  same  time 
provides  rather  simple  expressions  for  the  transition  rates. 

A  technique  developed  from  first  principles  and  not  restricted  to 

0 

perturbation  theory  has  been  published  by  Reiss. 

The  specific  process  studied  in  this  report  is  the  detachment  of 
an  electron  from  a  negative  ion  by  multiphoton  absorption.  It  will  be 
assumed  that  the  photon  energy  •h(u  is  significantly  less  than  the  detach¬ 
ment  potential  I.  Our  method  i.  not  applicable  for  evaluation  of  ioni¬ 
zation  of  neutral  atoms  or  molecules  unless  the  Coulomb  potential  in 

2 

the  final  state  is  screened  at  distances  less  than  (tU’j/c  ).  Study  of 
the  failure  of  our  theory  for  ionization  processes  demonstrates  what 
appear  to  be  significant  inaccuracies  in  the  theories  of  both  Kelaysh 
and  Reiss. 

Certain  approximations  will  be  made.  These  are: 

(1)  The  electromagnetic  field  of  the  laser  is  treated  classically 
and  is  assumed  to  vary  sinusoidally  with  the  frequency  w. 

This  implies  that  the  beam  is  intense  enough  that  quantum 
fluctuations  of  th-  E  and  B  fields  are  negligible  and  also 
that  the-  beam  coh-jrance  time  is  long  compared  with  that  for 
the  transition  to  occur. 

(2)  It  is  assumed  that  only  electric  dipole  transitions  need  be 
considered,  that  the  phQton  wavelength  is  very  large  compared 
wit’i  atomic  dimensions,  and  that  relativistic  corrections  for 
electronic  motion  may  be  neglected. 

(3)  We  shall  assume  that  only  a  single  electron  orbital  is  involved 
in  the  transition.  That  is,  we  shall  treat  the  problem  as 
though  it  were  the  case  of  a  single  electron  moving  in  a 
specified  binding  potential.  We  thus  avoid  the  complexity 

of  many -body  interactions. 


II  FORMULATION 


Let  h  =  h(p,x)  represent  the  Hamiltonian  of  the  atom  (or  molecule). 
Here  we  suppress  the  explicit  dependence  on  all  coordinates  except  the 
momentum  p  and  space  coordinate  x  of  the  active  electron.  In  the 
presence  of  the  electromagnetic  field  the  full  Hamiltonian  is 

H  =  h(p  -  a,  x)  CD 


where 


— *  e  -» 
a  =  —  A 
c 


(2) 


with  A  the  vector  potential  of  the  electromagnetic  field.  We  shall 
treat  the  vector  potential  as  classical,  writing  it  in  the  form 


cos  (cut  +  0  ) 
o 


(3) 


where  <u  is  the  angular  frequency  of  che  field  and  0^  is  a  constant. 
We  note  that 


-a  =  eE  =  eE  sin(uut  +  0  ) 

o  o 


/  A\ 

eE  =  u>  a  <4> 


Here  E  Is  the  electric-field  intensity. 


The  Hamiltonian,  Eq .  (1),  may  be  written  as 


with 


H  =  H  +  A 
o 


(5) 


(6) 


The  quantity  H  represents  the  "unoerturbed'’  Hamiltonian,  and  A  the  part 
o 

that  leads  to  transitions.  It  would  evidently  not  be  sensible  to  treat 
the  term  a2/2m  as  an  interaction,  since  this  C-number  commutes  with  the 
electron  variables. 

The  initial  state  of  the  ion  (or  atom)  with  the  active  electron 

bound  to  it  is  0  ,  where 
a 


h  0  =  W  0 
a  a  a 


The  solution  of  the  time-dependent  Schrodinger  equation 


H  0  =  i  0 
o 


corresponding  to  the  initial  state,  is  then 


x  -  -iW  t  -iy(t) 
0(t)  =  0  e  a  e 
a 


where 


Y(t)  a  f a(t,)2/2mdt ' 


(7) 


(8) 


(9) 


(10) 


The  final  state,  when  the  electron  is  out  of  range  of  the  potential  of 

f 

the  residual  ion  (or  atom),  is  0  .  This  satisfies  the  Schrodinger 
equation 

[(p  -  a)2/2m]  0f  =  i0f  ,  (11) 

or 

f  -3/2  ik*x  -iv  (t)  (12) 

0  =  (2  nr)  e  e  k 


where 


Y.(t) 

k 


(13) 


The  complete  process  is  described  by  the  wave  function  ty(t)  that  satis¬ 
fies  the  full  Schrodinger  equation 


Ht  =  i* 


(14) 


with  the  boundary  condition  that 


Aim  <!f(t)  =  0f(t)  .  (15) 

t“>  ♦  ^ 


If  we  imagine  the  electromagnetic  field  to  be  turned  on  very  slowly 

in  the  remote  past  and  off  very  slowly  in  the  remote  future,  the  S-matrix 

* 

for  the  transition  is 

S  =  Aim  ^f(t),  0(t)^  .  (16) 

t-*+a> 


See,  for  example,  M.  L.  Goldberger  and  K.  M.  Watson,7  Note  that  the 
outgoing-plane-wave  boundary  condition  is  used  in  Eq .  (16). 
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Following  standard  wave-packet  arguments,  we  have 


00 


—00 


A  perhaps  more  useful  form  for  the  S-matrix  is  obtained  on 
the  gauge  transformation 


*  = 


-»  -> 
ix«a  , 
e  ;• 


This  transforms  Eq,  (14)  into  the  form 


(h  +  F)jr  =  if 


where 


F(t)  =  -ex  .  E(t) 


(r 


making 


(11 


(1! 


(2 


The  S-matrix  in  this  representation  is 


The  total  rate  of  ionization  is  obtained  from  the  expression 


I  ~  2 


Here  T  is  some  long  time  interval, 
independent  of  T. 

The  expressions  (21)  and  (24) 
cussion  of  multiphoton  ionization. 


2 

dt((i,  A0)  (23) 

2 

dt(Y,  F$)  .  (24) 

chosen  large  enough  that  R  is 

were  used  by  Keldysh6  in  his  dis- 
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Ill  THE  QUASI STATIC  APPROXIMATION 


In  this  section  we  discuss  a  low-frequency  approximation,  assuming 

that 


■huu  «  I  =  -W 
a  a 


(25 


The  condition  (25)  was  assumed  also  by  Reissr  and  by  Keidysh.6 

The  limit  (25)  suggests  that  we  begin  with  solutions  to  the  eigen¬ 
value  equation 


where 


consistent  with  Eq . 


(15). 


(2n) 


-3/2 


ik-x 

e 


Jlimx-*» 


The  transformation 


(26 


(27 


(28 


converts  Eq.  (26)  into  the  form 


The  boundary  condition  (27)  leads  us  to  conclude  that 


W 

k 


=  K  /2m  s  e. 


so  the  0  satisfy  the  Lippmann-Schwinger  equation 

•4 

K 


+ 


(29) 


(30) 


If  a  can  be  treated  as  slowly  varying,  one  might  then  expect  that 


e.-  X 

-» 

k 


-iY,  (t) 
e  k 


(31) 


e 


-iY,  (t) 
k 


(32) 


would  provide  approxim>  <e  solutions  to  the  time-dependent  Schrodinger 
equation.  The  conditxons  under  which  Eqs.  (31)  and  (32)  are  valid 
approximations  will  be  studied  in  the  next  section. 


If,  for  the  moment,  we  assume  that  Eqs.  (31)  and  (32)  are  valid, 
the  expression  (21)  takes  a  rather  simple  form.  We  have 


(33) 
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Now 


-*  =  (0  ,  ex0  \ 

d  U  v 


is  just  the  dipole  matrix  element  for  the  single-photon  photoelectric 
process.  Thus,  knowledge  of  this  matrix  element  will  enable  one  to 
evaluate  the  multiphoton  transition  rate  when  Eq .  (32)  is  valid. 

For  evaluation  of  the  expression  (33)  the  Fourier  series  expansion 
of  Keldysh  is  conven  ent .  First  we  define 


2n 

U> 


3  =  — —  /  dt  sin(uut  +  0  )exp  ji[l  r  +  Y  (t)]| 
o  2«  /  o  L  a  k  J 


X 


/ 0  ,  ex0  \  •  E 

I  -*  al  o 

\  K  / 


(34) 


to  obtain 


2  2 
e  E 

S-  +2ni  Y  6  (i  +  e  +  - °  -  v-h  uj|S 

a— *  \  a  k  _ 

v=-oo  \  4mo) 


(35) 


Equations  (34)  and  (35)  were  used  by  Keldysh,  but  with  0_^  replaced 

K 

by  a  plane  wave: 


-  -3/2  it-x 

0  »  (2  tt)  e 

K 


(36) 
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IV  STUDY  OF  THE  QUASISTATIC  APPROXIMATION 


In  this  section  we  study  some  conditions  under  which  the  quasi^tatic 
expressions  (31)  and  (32)  might  be  expected  to  be  valid.  For  this  pur¬ 
pose  we  write 


=  X 


-iY,  (t) 

e  k 


(37) 


and  require  that  this  be  an  exact  solution  of  Eq .  (14).  On  noting  the 

relations  (4),  we  see  that  the  general  form  of  X  is 

k 


ix •  a 

X  =  e  0  (x,  a,  E)  .  (38) 

k  K 


Indeed,  substitution  into  Eq.  (14)  gives  us  the  following  equation  for 


0  r 

I? 


0 

K 


0 

•4 

K 


-  eE’V  0 
-»  -> 
a  K 


(39) 


First,  let  us  consider  the  strong  field  limit  and  assume  that  a 


can  be  treated  as  large.  Then  we  expect  the  Born  approx. mat  ion  to  be 

2 

valid  with  h  replaced  by  p  /2m. 
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In  this  case 


0 

K 


. -3/2  iK*x 
( 2  rr )  e 


(40) 


is  an  exact  solution  of  the  modifitd  Eq .  (39)  and  the  adiabatic  approxi¬ 
mation  is  valid.  Equation  (40),  substituted  into  Eq .  (34),  leads  just 
to  the  expression  used  by  Keldysh. 

The  above  argument  is  not  rigorous,  since  a  is  'arying  sinusoidally 
and  must  pass  through  zero.  It  does  make  it  plausible,  however,  that 
Eq.  (40)  can  represent  a  fair  approximation  when  (el  /ta)  is  large.  This 

O 

is  consistent  with  the  observation  of  Keldysh  that  .n  this  limit  (and 
for  u)  -»  0)  Eqs.  (40)  and  (35)  lead  to  Oppenheimer  ’  s  expression  (to  within 
a  numerical  factor)  for  ionization  of  a  neutral  atom  4v  a  static  electric 
field. 

We  next  use  Eq .  (39)  to  study  the  adiabatic  approximation  in  the 
weak  field  limit— that  is, 

|  a  |  «  k 
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Since  e  «  I  ,  we  take  the  eikonal  trajectory  to  be  a  straight  line 
k  a 

falling  into  the  force  center.  This  we  take  to  be  a  static  potential 
[  -V(x) ].  The  quanticy  Q  in  Eq .  (4?)  is  taken  to  be  the  positive 
root  of 


Q2  =  K2  +  2mV  +  eE*L^ (x)  +  a*L2(x)  (43) 

where  ~L  and  L  are  functions  to  be  determined.  To  ensure  convergence 

12  -1 
of  Eq.  (42),  it  is  assumed  that  V  (x)  decreases  faster  than  x  as  x-*». 

Thus,  if  V  is  a  Coulomb  potential,  we  assume  it  to  be  screened  at  large 

d.'Jtances.  The  quasistatic  approximation  is  valid  when  and  may  be 

neglected  in  Eq.  (43).  When  a  =  0,  Q  becomes 

r  2  -» i 1/2 

Q  =  j^k  +  2mV  (x)J  .  (44) 


If  we  let  the  eikonal  trajectory  lie  along  the  Z-axis  of  a  rectangular 
coordinate  system,  then  Q  =  Q(Z)  and  Q  =  Qq  (Z)  are  functions  of  Z  only. 
Then 

Q  (Z)  =  kjl  +  V  (Z)/V  (ZQ>]1/2  (45) 


where 


V(Z  )  = 
o 


(46) 
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io  first  order  in  a  and  E  we  obtain 


r  =  r 


o 


-> 

-  a 


a  *x 


L  dz7(2Q  ) 
1  o 


L  dz7(2Q  )  -  P 
2  o  a 


(47) 


where  F  is  the  function  T  for  a  =  E  =  0  (hence,  K  =  k)  and 
o 


(48) 


Substitution  of  Eq .  (42)  into  Eq .  (39),  with  the  standard  WKB 
approximation  of  neglecting  gradients  of  Q,  provides  us  with  the  coupled 
equations  (valid  to  first  order  in  a  and  ?) 


L  /(2m)  =  P  - 
1  a 


/L  f (2Q  )  dz' 
2  ° 

Z 


.  2 

L  /2m  =  u) 
2 


/L  (2Q  )  dZ# 
.  ° 


(49) 


where 
Eq.  (47) 


Use  of  these  relations  lets  us  re-express 


r  =  r 


a*x  +  a*L./(2m)  -  elf’L  /(2mo>  ) 


(50) 


The  quasistatic  approximation  is  seen  to  bo  valid  when 


L  S'  0,  L  /(2m)  S  P  .  (51) 

2  1  a 


To  integrate  Eq .  (49)  we  introduce  the  new  variable 


v(Z)  2 


/ 

o 


(k/Q  )dz7a 
o  o 


(52) 


Then 


_d_ 

dv 


-ofiL1/(2m) 


_d_ 

dv 


dP 

— -  +  a  L  /(2k) 
d  v  o  2 


(53) 


where 


fi  =  mcua  /k  =  t.'Jo/ 

o 


l2 


(54: 


The  quantity  (48)  may  be  expressed  in  terms  of  v.  We  first  define 


P(Z)  =  a  v(Z)  -Z 

o 


Then 


P  (Z)  =  P(®)  -  P(Z) 
a 


Equation  (53)  may  be  integrated  in  the  form 


L  /(2m)  =  P  -  H 
1  a 


/ 


sin  [fKv'  -  v)]P  (v/)dv/ 
a 


'55 


(56 


w 

L  /(2m)  =  (tu/n)  f  I1  -  cos[n(v'  -  v)]|r(v/)dv/ 


(5' 


Here 
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L  (Z)/(2m)  =  Li(0)/(2m)  +  Z  +  .  .  . 

L  (Z)/(2mu))  =  L  (0)/(2mu)> 

2  2 

-  -  (Z/a  )3/2  (a  /Z  )1,  2  n  L  (0)/ (2m) 

3d  o  o  1 

+  .  .  .  .  (61) 


The  linear  Z-term  in  L  corresponds  to  the  Reiss  factor  exp(i  a'x). 

The  L  (0)  and  L  (0)  terms  (independent  of  the  space  coordinate)  arise 
1  2  ^ 

in  satisfying  the  boundary  condition  at  |x|  =  ®.  These  are  time-dependent, 
so  they  will  contribute  to  Sq  [Eq .  (34)].  Indeed,  since  1^(0) /2m  »  aQ, 
we  see  no  reason  to  neglect  these  terms. 

When  the  screening  radius  is  R  «  Z^  (but  »  a^)>  we  can  evaluate 
the  relative  error  in  the  quasistatic  approximation  from 


r  =  L 

’i  L  i 


(0)/ (2m)  -  P&(0)] 


/P  (0) 
a 


(R/a  )3(hu)/Ry)2/(12)2  .  (62) 

o 


When  this  quantity  is  sufficiently  small  compared  to  unity,  we  expect 
the  quasistatic  approximation  to  be  valid. 

For  R  »  Z  the  quasistatic  approximation  does  not  appear  to  be 
o 

valid.  Numerical  evaluation  of  Eq .  (62)  gave  values  of  gj^  and 

g  =  |l  (0)/[2mu>  P  (0)]|  (63) 

2(2  a  i 


to  be  about  unity  for  a  wide  range  of  parameters.  Furthermore,  an 
estimate  of  the  second-order  terms  in  a  and  ?,  omitted  from  Eq.  (50), 
indicated  a  contribution  comparable  to  the  first-order  terms  when  Sq 
[Eq .  (34)]  is  evaluated. 
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We  conclude ,  then,  that  for  a  final  ionic  state,  with  rcreening 

distance  R  »  Z  ,  we  are  not  able  to  evaluate  multiphoton  ionization 
o 

rates.  Since  we  have  found  what  appear  to  be  large  corrections  in  the 
expressions  used  by  Keldysh5  and  Reiss,6  our  confidence  is  not  great 
in  the  accuracy  of  their  results. 

We  are  able  to  draw  a  more  optimistic  conclusion  for  detachment 
of  an  electron  from  a  negative  ion  when  the  residual  molecule  does  not 
have  p.  significant  permanent  electric  dipole  moment.  In  this  ase  the 
potential  in  the  final  state  has  the  asymptotic  form 

-V(r)  =  -  a^e  /  (2r  ) 

r  -  *  co 

where  a  is  the  polarizability. 

P 

-  3 

For  photodetachment  of  0  \vr  have  a  w  5.29  and  1^  =  1.46  eV, 
Equation  (67)  should  be  numerically  evaluated  for  a  potential  with  the 
asymptotic  form  (64)  and  for  small  r  modified  to  remove  the  singularity 
at  r  =  0.  The  error  functions  g  and  g  are  shown  in  Figure  1.  The 

i.  & 

curves  are  labeled  by  N,  the  least  number  of  photons  required  for 

detachment,  and  are  given  as  functions  of  e  =  N  liu'  -  I  .  As  anti- 

cipated,  we  see  that  for  large  N,  and  an  not  too  small,  both  g^  and 

g  become  small,  as  required  for  the  validity  of  the  quasistatic 
2 

approximation. 
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drawing. 


V  DETACHMENT  OF  0 


We  now  use  the  expression  (34)  to  evaluate  the  m»ltiphofon  detach¬ 
ment  rate  of  0  in  the  weak  field  limit.  The  single  photon  detachment 
rate  has  been  studied  experimentally  by  Branscomb  et  al.9  Theoretica’ 
studies  have  been  made  by  Brueckner  and  Klein10  and  by  Gillespie.11 
Since  the  attached  electron  is  in  a  p-orbital  state  the  transition  will 
take  place  to  a  state  of  zero  angular  momentum  if  the  electron  energy 
is  sufficiently  low.  We  thus  use  condition  (25)  to  write11 


f_1(-K)  CE 

o 


(65) 


where  C  is  a  constant  and  f  is  the  Jost  funrtio'  ,hich  we  write  in  the 

♦ 

effective-range  approximation  as 


f 


-1 


(-K)  =  (K  +  ia)/(K  +  ib) 


(66) 


Here  a  and  b  are  constants  related  to  the  scattering  length  and  effective 
range  for  electron  scattering  by  a  neutral  cxygen  atom.^ 


See,  for  example,  Ref.  7,  p.  540. 

T.  F.  O'Malley  et  al.12  have  shown  that  for  a  potential  having  the 
asymptotic  form  (64)  the  effective  range  expansion  is  modified.  For 
the  very  low  electron  energies  of  interest  to  us  this  modification 
is  not  significant. 
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Following  Keldysh,5  we  take  0  =  —  and  introduce  the  variable 

o  2 

u  =  sin  ast  in  Eq .  (34).  We  shall  assume  that  the  electromagnetic  field 
is  weak  in  the  sense  that 


(eE  /(u)  «  k 
o 


(67) 


We  rewrite  (34)  in  the  form 

■„  -  -  «*>-1  Lfe  eiQ<k'u>]  °Eo  • 

Here  the  clos  contour  encloses  the  points  u  =  -1,  1  and  lies  within 

iQ 

any  other  singularities  of  e  /f . 

In  the  weak  field  limit  we  can  deform  the  contour  so  that  every¬ 
where  on  it  |u|  »  1.  An  elementary  evaluation  gives  us,  then, 

iQ  N  2  r  -»  -»  2  t  — N 

e  =  (i/2)  exp(-  0u  /4)  exp  |-ek -E^/ (mu)  )Ju  ,  |u|  »  1  .  (69) 


Here 


N  = 


2  2  2  "" 
e  E  /  ( 4  ukd  ) 


/(hu)) 


is,  according  to  Eq.  (35),  the  number  of  photons  required  to  detach 

[  2  2  2  1 
e  E  /(muj  )J  /Chou). 

As  Keldysh  observed,5  Eq .  (69)  has  two  saddle  poiifrts  at  the  roots 
of 


I 

a 


0 


(70) 
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These  correspond  to 


u 

s 


-i-y  = 


*iu> 


/(eE  ) 
o 


(71) 


As  also  observed  by  Keldysh,  the  expression  (65)  has  poles  at  the 
positions  defined  by  (70).  The  quantity  C.  in  Eq.  (65),  which  is 
really  K-dependent,  has  poles  at  the  same  positions.11 

Keldysh  used  the  Born  approximation,  so  in  his  expressions,  f  =  1. 
He  was  thus  able  to  evaluate  (68)  keeping  only  the  contribution  from 
the  saddle  points  at  u^  Sr  tiV. 

We  cannot  use  the  Keldysh  approximation  because  f  (-K)  has  branch 
points  at 

u  =  uj(-k||  ±ik^)/(eE^)  <72 


where  k)(  is  the  component  of  k  parallel  to  and  kA  is  the  component 
of  k  perpendicular  to  Eq.  Since  |uj  «  |uj  our  integral  (68)  will 
receive  its  principal  contribution  from  the  vicinity  of  the  branch 
points.  The  contour  of  integration  used  is  shown  in  Figure  2. 

Evaluation  of  (68)  is  straightforward  when  N  is  large.  The  rate 
of  detachment  is,  finally,  in  the  limit  of  large  N: 


=  2n  1 

r  3  ( 

-t  k6  (i  + 

n  2  2 

e  +  e“E  / (4muj  ) 

-  Nhu))  |  £ 

J 

\  a 

k  o 

r  ,  v  1 2 

/ 

r  “l 

e. 

(1) 

2" 

a, (k)/a 

L  1  o  . 

1  L  1  l\\ 

k 

Sh  — 

'  2jt 

|  L  \b  “  I/J 

3 

2oN  Ry_ 

U(2e^/hou)  |eEQ/(2ku))^ 


2N 


(73) 


*  Equation  (71)  would  be  exact  if  £  were  zero.  For  nonvanishing  k  the 
actual  saddle  points  are  displaced  slightly  from  the  position  ±  i  Y. 
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FIGURE  2  CONTOUR  OF  INTEGRATION  FOR  THE 
INTEGRAL  (68),  SHOWING  THE  TWO 
BRANCH  POINTS  ub  AND  THE  TWO 
SADDLE  POINTS  us.  Branch  lines  are 
indicated  between  u  =  -1  and  +1  and 
from  ±ub  to  ±us. 

In  tie  lower  version  of  this  equation,  k  is  evaluated  when  the 

argument  of  the  6-function  vanishes,  or  is  the  fine-structure  constant 

and  R  is  the  Rydberg  constant.  We  have  expressed  R  in  terms  of  the 

cross  section  a(k)  for  detachment  by  a  single  photon  of  energy 

■ho)  =  I  +  e  .  Finally, 

1  a  k 


It  is  interesting  to  observe  that  for  large  photon  numbers  the 


cross  section  (73)  is  considerably  enhanced 
Our  weak  field  approximation  (67)  fails  for 
however,  and  the  apparent  divergence  of  Eq . 
failure  of  this  approximation. 


for  electrons  near  threshold, 
very  low  electron  energies, 
(73)  at  k  =  0  is  due  to  the 
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VI  DETACHMENT  OF  0 

2 

Cross  sections  for  photodetachment  of  0^  have  been  published  by 

Burch  et  al.13  The  detachment  energy  is  I  =  0.44  eV  and  the  electronic 

a 

2 

configuration  is  a  II  state.  An  electric^dipole  transition  cannot 

e  2 

take  place  to  a  final  S-state  for  the  II  state.  Geltmann14  has  dis- 

g 

cussed  the  implications  of  this  for  the  near-threshold  dependence  of 
the  single-photon  detachment  cross  section. 

A  theory  of  final-state  enhancement  does  not  seem  to  have  been 
given  for  nonspherical  potentials.  Until  such  a  theory  is  available, 
we  shall  use  Eq .  (73)  with  an  appropriate  shift  of  parameters.  For 
oi<k)  we  shall  take  the  experimental  value  of  Burch  et  al.13  This  is 
admittedly  not  a  rigorous  theory.  If,  however,  the  dipole  matrix 
element  has  the  branch  points  (72)  as  closest  to  the  origin  in  the 
4-plane,  the  ionization  rate  is  rather  insensitive  to  other  details 
of  the  matrix  element. 

The  resulting  ionization  rate  as  a  function  of  laser  power  is 
shown  in  Figure  4  for  10.6-4m  radiation. 
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DETACHMENT  HATE 


FIGURE  4  THE  DETACHMENT  RATE  FOR  Oj  SHOWN 
AS  A  FUNCTION  OF  LASER-BEAM  POWER 
FOR  WAVELENGTHS  OF  10.6  AND  5.5  fim 
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